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I. INTRODUCTION 

In Refs. GJ it was demonstrated how a quantum com- 
puter could be built using only single photon sources, pas- 
sive linear optics elements, and photon detectors. Quantum 
computing proposals that use photons to encode information 
are particularly interesting because of practical applications 
to quantum communication over optical fibers, and the natu- 
ral resilience of photons to decoherence. This proposal is also 
a conceptual departure from other quantum computing pro- 
posals because it requires post-selection of states in order to 
overcome the limitations imposed by the choice of physical 
resources - namely, the fact that one cannot make photons in- 
teract using passive linear optics elements. Moreover, even 
when considering ideal hardware (i.e. lossless and infinitely 
precise linear optics elements, 100% efficient detectors) one 
must use error correction codes to make the implementation 
efficient 01 ■ However, photon detectors are necessary for the 
qubit measurements used in post-selection and error correc- 
tion, and good photon detectors are notoriously hard to build. 
In this paper, we investigate the maximum error rate that a 
linear optics quantum computer, as proposed in Ref. yj can 
sustain, assuming that the only source of hardware imperfec- 
tions is the finite photon loss at the photon detectors. 

The paper is organized as follows. In Section|ll]we briefly 
describe the construction of a probabilistic gate as given in 
Ref.Q] and emphasize the description of the error model which 
arises naturally from considering those gates. In Section Hill 
we describe an error model based on the same gate construc- 
tions, but assuming that the single photon detectors have less 
than perfect efficiency. In Section HVl the error correction code 
and the circuits used for error correction are described, along 
with constraints for fault-tolerance. Finally, in Section lvTll the 
recursion relations for the error rates at different level of con- 
catenated encoding are given under worst case assumptions, 
along with a brief description of how they were calculated, 
and the threshold values are stated. 

Throughout the paper we will use the Pauli matrices X = 



( i o )' Y = ( i o )' z = ( o -? ) along with the 

identity matrix 1 . 



II. IDEAL HARDWARE ERROR MODEL 

In the efficient linear optics quantum computing proposal 
put forward by Knill, Laflamme and Milburn |1], qubits are 
encoded as a single photon in one of two optical modes, that 
is, |o) is represented by the photon number state |0) and 
|i) is represented by the photon number state |1) |0) - this en- 
coding is also called the dual-rail encoding 12611 . The only re- 
sources available are single photon sources, passive linear op- 
tics elements (such as beam splitters and phase shifters), and 
photon detectors. In our model, the only source of hardware 
imperfection is the efficiency of the photon detectors, that is, 
the photon sources and passive linear optics elements are as- 
sumed to be perfect. We also assume that classical computa- 
tion and control are delay and error free, and that all sources 
of failure - teleportation and measurement failures - are sta- 
tistically independent. 

While single qubit operations can be efficiently performed 
using only phase shifters and beam-splitters {51 - and there- 
fore can be considered error free in our model - two qubit 
operations require state post-selection through measurement 
of ancillary modes. If the desired measurement is not ob- 
tained, the operation may or may not have been applied prop- 
erly, but if the desired measurement is obtained, the proper 
operation is guaranteed to have been applied. In Ref. Q] a 
probabilistic sign shift gate, which performs the operation 
NS_ = |0}(0| + |1)(1| - |2}{2| on number states of a given 
mode, is described. This gate succeeds with probability j. A 
probabilistic entangling gate CSIGN, which performs the op- 
eration |oo) (oo| + |oi) (oi| + |io) (io| — |n) (n| on two qubits, 
can be constructed by using two NS_ gates along with two 
beam splitters, but no extra ancillae or measurements. Since 
both NS_ applications must succeed, the overall probability of 
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FIG. 1: Gate teleportation protocol for a CSIGN. The resource state 
is | Boo) = and the two qubit measurement is a Bell basis 

measurement. The double lines indicate operations conditioned on 
the outcome of the measurement. 

success for a CS IGN is jg . In order to make this construction 
scalable, one has to use gate teleportation |4|. which turns 
the gate construction problem into a state preparation prob- 
lem with the advantage that such states can be prepared off 
line without a priori knowledge of the inputs to the gate. The 
idea of gate teleportation is to use the conjugation relations of 
certain gates to modify the resource state and the correction 
operations of the teleportation protocol, so that these gates 
can be applied implicitly during the teleportation, in a manner 
similar to the protocol for the CSIGN gate illustrated in Fig- 
ure ^ Strictly speaking, the CSIGN construction is based on 
the teleportation of some of the modes constituting the qubits 
that must interact, and thus the operations are not general uni- 
taries but instead restricted to linear optics operations, single 
photon sources, and post-selection based on single photon de- 
tection of certain modes. 

Even though the direct C S I GN construction has a small (but 
non-zero) probability of success, it is always known when the 
teleportation succeeds, and the resource state of the teleporta- 
tion is independent of the inputs to the operation so one can 
make many attempts to produce such a resource state before 
performing the teleportation. Alternatively, one may make 
many attempts in parallel, and select one that succeeds. While 
without teleportation the number of attempts needed for er- 
ror free quantum computation grows exponentially for circuits 
constructed with the probabilistic CSIGN gates, with the gate 
teleportation it grows linearly QQ]. Since the CSIGN along 
with all single qubit unitaries form a universal set for quantum 
computation |5], this gives an efficient construction for a uni- 
versal quantum computer using only linear optics elements. 

Teleportation in linear optics cannot be performed deter- 
ministically 0], but it can be performed with probability ar- 
bitrarily close to one at the cost of higher complexity for the 
teleportation protocol 01- Once again, as in the probabilistic 
CSIGN construction, one can determine whether the proto- 
col succeeded through the measurement of ancillary modes, 
so possible failures are always flagged by the measurement 
outcomes. Since it is not clear how to increase the probabil- 
ity of success of the direct CSIGN construction, and bounds 
for possible probabilities of success are predicted to be signif- 
icantly smaller than one |7J], this gate teleportation scheme is 



highly necessary. We say that the probability of these intrin- 
sic teleportation failures occurring is e;d e ai, since such failures 
occur even when considering ideal photon detectors (i.e. with 
perfect efficiency) and linear optics elements. 

What makes the teleportation scheme highly attractive as 
well is exactly what the failure entails. When one of the tele- 
portations involved in the CSIGN fails (but no detectors fail), 
the output is equivalent to a successful teleportation followed 
by a measurement in the Z eigenbasis - that is, we can view it 
as a successful gate application followed by a Z measurement 
of one of the qubits independently with probability eideai, and 
such an event is automatically flagged by the outcome of the 
measurement of the ancilla modes. This is what we call the 
ideal hardware error model, and it has been shown to be very 
benign, with an error threshold arbitrarily close to one 0IH 

Hence teleportation failures under this model mean that 
one of the projectors into the eigenbasis of Z has been ap- 
plied. Exactly which projector was applied to the teleported 
qubit depends on the outcome of the measurement of the an- 
cillae. Formally, the projectors are Z + = i (1+Z) and 
Z = j (1 — Z). Since the location of these measurements 
is flagged by post-selection, this type of failure is a form of 
erasure - an error of known location. 



HI. LOSSY HARDWARE ERROR MODEL 

One of the largest technical hurdles in the implementation 
of these probabilistic gates is the fact that single photon detec- 
tors are notoriously difficult to build. While for some wave- 
length very high efficiency can be obtained, the rate at which 
false photon detections are signaled, the so called dark counts, 
is unacceptably high \ 9]. Dark counts are particularly trouble- 
some for the gates proposed in Ref.Q]because they could cause 
photon loss to go undetected, as well as causing incorrect 
post-selection of states. Recent proposals of photon detectors 
based on phase transitions in superconductors have very low 
dark count rates, although little attention has been given to 
optimizing the efficiency of these detectors. We are interested 
in finding the minimum efficiency necessary for these detec- 
tors, given dark counts stay negligible, in order to be able to 
perform useful quantum computation 12711 . 

Modifications to the gate teleportation protocol that allow 
for the detection of photon loss at the detectors are known (]]]. 
This protocol can differentiate between the teleportation fail- 
ures due to the limitations of linear optics (leaving the error 
model due to such failures intact) and the failures due to pho- 
ton loss at the detectors - when this occurs, the corrupted qubit 
is replaced with a fresh qubit in a known state. We take the de- 
tection of a single physical qubit to fail with probability 6, and 
the overall probability of the CS I GN gate teleportation failing 
due to photon loss at the detectors to be ei oss . The form of 
ei oss as a function of 6 depends on the choice of protocol used 
- Ref. Q] describes a family of protocols - and for the purposes 
of this paper, we will take them to be independent parameters. 
Once again, any possible error due to these types of failure is 
always flagged by the gate construction, which significantly 
simplifies error correction since it is known a priori where the 
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failures have occurred. 

Consider the error model due only to photon loss at the 
detectors (ejdeai = 0) - such an error model is not physical, 
since a linear optics quantum computer will always have tele- 
portation failures, but taking this limiting case simplifies the 
analysis significantly. If we consider a single physical qubit 
measurement where photon loss occurred, it is clear that all 
information about the qubit is lost. We can model this loss of 
information by the full erasure superoperator 

E(p) = \t = J(p + XpX + ZpZ + YpY). (1) 

However, because of the dual-rail encoding, it is always clear 
which physical qubit measurement failed in such a manner 
since at the lowest level of encoding qubit measurement con- 
sist of measuring the two constituent modes, and thus only 
one of the detectors may click. If a photon is not detected on 
either mode, a measurement failure has occurred, and we may 
replace the qubit with a fresh qubit in a known state. This 
is fundamentally different from the depolarizing channel Jioll 
where there is no a priori knowledge of the position of the 
errors - one may think of this superoperator as the depolariz- 
ing channel superoperator conditioned on perfect information 
about which qubits were randomized. 

It is clear that the full erasure does not commute with gate 
teleportation, otherwise we could in principle transfer some of 
the information from the control qubit of a CSIGN to the tar- 
get even though all information was lost in the failed telepor- 
tation of the control qubit. We take a worst case approach, and 
assume that any photon loss during one of the teleportations 
yields total information loss of the qubit being teleported, and 
that the error model for the other qubit associated with the 
same gate teleportation in which the photon was lost is de- 
termined by which correction might have been needed to be 
applied. In the particular case of the linear optics proposal, the 
error model is symmetric: if we disregard the classical corre- 
lation between the errors of the outputs of the CS IGN, photon 
loss in one of the teleportations translates to a full erasure of 
the qubit being teleported and a Z erasure of the other qubit 
involved in the gate operation. This is because in both cases 
the correction operation from one teleportation to the other is 
a Z gate, and if it is not known whether such gate was sup- 
posed to be applied because of the photon loss, we have the 
superoperator 

Z(p)=-(p + ZpZ), (2) 

which is what we call a Z erasure. Note that this superoper- 
ator can be interpreted as an unintentional Z measurement of 
unknown outcome, since 

Z +P Z++Z_pZ_ = l(p + ZpZ) =Z(p), (3) 

and thus this is fundamentally different from the phase erasure 
channel fTofc . since it provides the added a priori knowledge 
of when corruption has occurred or not. As mentioned be- 
fore, we say that there is a probability ei oss that photon loss 



occurs in one of the teleportations in the C S I GN implementa- 
tion, which entails a full erasure of the qubit being teleported, 
and a Z erasure of the qubit to which it was coupled through 
the CSIGN application. The correlation between the erasures 
on these two qubits is ignored for our calculations, but could 
be exploited to obtain better thresholds. 

In summary, all single qubit operations are taken to be error 
free, and the CS IGN is taken to introduce at each output qubit 
either Z measurements, full erasures, or Z erasures with fi- 
nite probability. For simplicity, we consider the different error 
models independently, that is, we calculate the threshold for 
the case where only Z measurements occur (where the hard- 
ware is ideal, with ei oss = and 6=0, but teleportation 
is imperfect, i.e. eideai ^ 0), and we calculate the threshold 
where only full and Z erasures occur (that is, where the tele- 
portation protocol is perfect, £id ea i = 0, but the detectors are 
not, so 6 7^ and ei oss ^ 0). 

A collection of single qubit erasures is referred to as an era- 
sure pattern, and the weight of the erasure pattern is the num- 
ber of qubits that have been affected by an erasure, regardless 
of the type of erasure. 

IV. FAULT-TOLERANCE AND ERROR CORRECTION 

The quantum error correction code considered here to pro- 
tect the data from the error model in question is a [[7, 1 , 3]] 
self-orthogonal, doubly-even CSS code ifTTl fl2il with stabi- 
lizer generators 

Mi=X®X®X®X®1®1®11 

M 2 =X®X®1®1®X®X®11 

M 3 =X®1®X(X)1®X®1®X 

M 4 = Z(g)Z(g)Z®Z®ll®]l®l ( ' 

M 5 =Z®Z®1®11®Z®Z®]1 

M 6 = Z®l®Z®l®Z(g)l®Z. 

Self-orthogonal CSS codes are particularly suited for the er- 
ror model considered here because they yield simple construc- 
tions of fault-tolerant encoded Clifford group operations. Al- 
though the Clifford group is not a universal set of quantum 
operations, it is well known how to extend it in order to ob- 
tain a universal set fault-toleran tly Jl3l Il4ll . as well as how 
this affects the threshold value IU5I1 . In the case of the er- 
ror models considered here, the threshold is unaffected since 
we require that computation be performed only on error free 
states - that is, error correction is performed until an uncor- 
rectable error occurs, causing the computation to be aborted, 
or until the data is error free, at which point the computation 
may continue. 

Recall that the Clifford group consists of all operations that 
preserve the Pauli group under conjugation, and this group 
is generated by the CSIGN gate, the Hadamard gate H = 
^ ( l Ji ). and the phase gate P = (J °). The CSIGN can 
be implemented transversally by qubitwise CSIGNs between 
two encoded qubits. Note that CSS codes have a transver- 
sal encoded CSIGN. If in addition H has to be transver- 
sal, then the CSS code has to be constructed from a single, 
self-orthogonal classical code. Moreover, if P also has to 



be transversal, then all codewords in this classical code must 
have doubly-even weight. Since the 7-qubit code considered 
is both doubly-even and constructed from a self-orthogonal 
classical code, the Hadamard and phase gates can be im- 
plemented using only single qubit operations qubitwise, and 
therefore are error free under the models considered here. In 
order to obtain a universal gate set, we can add a non-Clifford 
gate such as the 7t/ 8 gate, which has a known fault-tolerant 
construction ffl . We will make use of the properties of 
doubly-even CSS codes for the error correction circuits stud- 
ied in Section lyfl 



V. ENCODED ERROR MODEL 
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FIG. 2: A circuit for correcting a single Z measurement in the top 
qubit as long as the remaining 3 qubits are erasure-free. The four 
qubits are assumed to be part of a block encoded with the 7 qubit 
code. The state |+) is the +1 eigenvalue eigenvector of X. The dou- 
ble lines indicate control based on the outcome of the measurement. 



Error models consisting of erasures yield particularly sim- 
ple encoded error models because it is always known when 
an error is unrecoverable. In such a case one can simply take 
the code block to be an encoded erasure at the next level of 
encoding. In general, an uncorrectable failure is not an en- 
coded failure, and it requires further processing to make it an 
encoded failure, as will be discussed below. What remains to 
be determined is what kind of encoded failure results. 

In the case of ideal hardware, all erasures are Z mea- 
surements of known outcomes. Turning to the 7 qubit code 
described above, we have that all weight one and two era- 
sure patterns are correctable. Of the 35 possible weight 
three erasure patterns, 28 are correctable. The remaining 7 
weight three patterns, along with all other patterns with higher 
weight I2H1 . can be identified with an encoded Z measure- 
ment, and are therefore uncorrectable failures. These collec- 
tions of individual single qubit measurements are not, strictly 
speaking, equivalent to an encoded measurement, since they 
collapse the seven qubits into a state outside the code space. 
However, because the encoded |o) and |i) are superpositions 
over mutually exclusive sets of states, it is easy to infer which 
encoded state the measurement results correspond to, and then 
replace the qubits with a fresh encoded |o) or |i). This opera- 
tion is taken to be error free, since we assume that state prepa- 
ration can be attempted until no errors have occurred. Since 
all encoded failures are Z measurements, in the ideal hard- 
ware error model the break even condition between encoded Z 
measurements at the first level of encoding, denoted by Z' 1 
and single qubit Z measurements, PrfZ' 1 ' measurement) = 
Pr(Z measurement), implies 

6 ideal = e ideal, (5) 

where ejjjj is the failure rate at the first level of encoding. 

In the case of the lossy error model, because each qubit 
can suffer either Z erasures or full erasures, different uncor- 
rectable errors on the 7 qubit code will lead to different en- 
coded errors. For simplicity we can take all encoded failures 
to be encoded full erasures, which in the first level of encod- 
ing we denote E' 1 ', so that the break even condition on the 
probabilities Pr(E ( 1 > ) = Pr(E) implies 

e ios! = 2 £loss ' ® 



since only half of the lossy error model failures are full era- 
sures, where ejJJ is the failure rate at the first level of en- 
coding. The resulting threshold is at most as high as the real 
threshold, considering the different kinds of encoded erasures 
that are simpler to correct, but should be lower in general. A 
more detailed analysis can be made Jl7lfl8ll . and exact prob- 
abilities distributions for the different kinds of erasures can 
be calculated using the same technique used to calculate the 
threshold in Section iyTD but the simplifying assumption made 
here is enough to match the prediction in Ref . y. 

In reality we would like to consider an error model that 
takes both these sources of error into account. Section IIVI 
describes how each of the different kinds of erasures are cor- 
rected, and demonstrates the progressively higher cost of cor- 
recting a Z measurement, a Z erasure and a full erasure. Given 
this fact, for the error correction code chosen here, and for the 
error correction technique used here, the threshold for an error 
model consisting of both types of failures is bounded above by 
the ideal error model threshold, and below by the lossy error 
model threshold. 



VI. ERROR CORRECTION CIRCUITS 

In general, in order to correct errors by using a stabilizer 
code, one simply needs to measure the stabilizer generators 
and infer the most probable error that occurred and apply the 
correction. In the case of erasure errors, the knowledge of 
which qubits have been affected by the error superoperator 
greatly reduces the number of stabilizer operators that need to 
be measured. This is because we need only measure stabilizer 
operators that act non-trivially on the qubits affected by the 
error superoperators, and this greatly reduces the probability 
of introducing more errors into the data. 

This procedure can be further optimized by incorporating 
the syndrome measurement and correction into a single step, 
as first described in Ref. for the case of Z measurement cor- 
rection codes. 

The fact that the 7 qubit code employed here is based on 
a classical doubly-even code, allows us to consider the 4 
qubit subsystem in the support of any given stabilizer operator. 
Since this is a CSS code, we can focus on stabilizer operators 
that are made up of tensor products of Xs and identities, and 
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FIG. 3: Performing a CNOT between |+) and an arbitrary state |ot) 
through teleportation. 



stabilizer operators that are made up of tensor products of Zs 
and identities. 

Since we are only interested in correcting erasures of 
weight up to 3, we can consider a single qubit that has un- 
dergone some erasure, along with three qubits that are still 
intact. We choose stabilizer operators that act non-trivially on 
all of these four qubits and trivially on all other qubits - this 
is always possible in the 7 qubit code. Considering only this 4 
qubit subsystem, the stabilizer operators in question are 1 29] 



M 1 = X ® X ® X 



M, 



Z <g) Z ® 



X 

Z tg> z. 



(7) 



The usual approach is to measure these two operators fault- 
tolerantly in order to determine what kind of Pauli correction 
needs to be applied to the erased qubit. If we consider Z era- 
sures, we need only measure Mj. Alternatively, one can sim- 
ply use the circuit depicted in Figure |2] (without loss of gen- 
erality, we consider the corrupted qubit to be the first qubit 
of the four) ll30ll . As demonstrated before, a code that can 
correct Z erasures can also correct Z measurements at known 
locations, since a Z erasure can be given by a Z measurement 
of unknown outcome. In this case, we do not need to perform 
the explicit measurement in the Z eigenbasis since we already 
have that information, but the rest of the circuit remains as in 
Figure|2]- thus there is an added cost of a possible measure- 
ment failure when correcting Z erasures. If the measurement 
does fail and the qubit is destroyed (as is the case when a pho- 
ton is detected), one simply abandons the attempt at correcting 
the Z erasure, since it will then be replaced by a full erasure. 

The circuit in Figure|2]is not fault-tolerant, as it stands. In 
order to make it fault-tolerant, one simply applies a modified 
teleportation protocol to the three error- free qubits (J, 0, fl7ll . 
In order to understand why this is fault-tolerant, consider the 
teleportation of the three intact qubits before applying the 
CNOTs in Figure |2 Since CNOT is a Clifford gate, we can 
apply the CNOTs to the Bell states needed for teleportation, 
and simply modify the recovery stage of the teleportation, in 
a manner similar to what was described in Figure^ Note that 
it is unnecessary to teleport the control qubit, i.e. the qubit 
in the state |+), since it is a fixed resource state. Unlike the 
C S I GN gate construction, where photon modes are teleported, 
this is a teleportation of the qubits, and can be thought of in 
terms of the usual higher level gates such as CNOTs and Pauli 




FIG. 4: Fault-tolerant measurement of stabilizer operator Z g> 1 (g> 
Z®ll®Z<g>ll®Zas proposed by Shor IT3l . The qubits in the 
shaded region are ancillas, and the measurement of the operator is 
inferred from the measurement of the 4 ancilla qubits. 



operators. 

This procedure becomes clearer if we consider the telepor- 
tation of only one of the qubits, followed by the CNOT with 
the resource state |+), and propagate the CNOT backwards, as 
described in Figure|5] If there is a failure in the Bell measure- 
ment of the qubits, it will only possible cause a Z erasure in 
the control bit, as well as a possible full erasure on the tar- 
get bit. Since Z errors do not propagate from the control of 
the CNOT gate, the other two CNOTs can be performed fault 
tolerantly in parallel. 

In summary, we simply need to measure the erased qubit 
in the Z eigenbasis, discard the erased bit, and apply the 
teleportation-based fault-tolerant version of Figure [2] to three 
erasure free qubits in the codeword and an extra qubit in the 
|+) state to replace the discarded one, and the Z erasure will 
be corrected. If any of the teleportations fail, the failure will 
affect only the qubit we are attempting to recover and the qubit 
that was being teleported. 

The circuit in Figure |4] is used to measure stabilizer opera- 
tors made of tensor products of Zs and Is, and thus it partially 
corrects full erasures, yielding a Z erasure if successful. The 
fault-tolerant version of the circuit in Figure[2]can be used to 
correct Z erasures as well as unintentional Z measurements. 
The strategy taken is to correct full erasures first, and once 
there are no more full erasures, to correct Z erasures and Z 
measurements. We do not claim that this strategy and circuits 
are optimal for error correction, but this choice simplifies the 
exact calculation of the threshold significantly. 



VII. THRESHOLDS 

The threshold theorem USUI EH states that by concate- 
nated coding - that is, the repeated encoding of a quantum 
state - one can perform quantum computation with arbitrarily 
small error efficiently as long as the error is below a certain 
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threshold. We take this threshold to be the smallest probabil- 
ity of error such that the probability of an encoded failure e' 1 ' 
is equal to the probability of a single unencoded qubit failure 

e'°> = e. 

The encoded error rate can be calculated by tracking the 
probability of going from any given erasure pattern to any 
other erasure pattern during an attempt at error correction. 
This describes a Markov chain, and such a description can be 
made more compact by considering symmetries of the error 
correction code and of the error correction circuitry fl7il . Era- 
sure patterns can be grouped into equivalence classes defined 
by the error correcting code as well as the error correcting 
procedure, and we need only consider probabilities of going 
from one equivalence class to another - in the case of the 7 
qubit code, we need to consider only 1 1 equivalence classes 
versus the 1 28 that would be necessary for a naive descrip- 
tion of the Markov chain. It is straightforward to obtain the 
initial distribution of the different equivalence classes as well 
as the transition matrix of the Markov chain that describes the 
change in the distribution due to one error correction attempt. 
The distribution after multiple error correction attempts can be 
obtained by taking higher powers of the transition matrix and 
applying it to the initial distribution. Each of the non-trivial 
equivalence classes of erasure patterns can be associated with 
an encoded erasure, and therefore one can obtain the error dis- 
tribution at any given encoding level. In some cases, there is 
additional processing and transitions between different equiv- 
alence classes associated with the mapping between the era- 
sure pattern and an encoded erasure, since, for example the 
erasure operator l<g>ll(g>ll<8>llCS>ll®ll£>5Eis not equivalent 
to any encoded operation. In order to account for such pro- 
cessing, another transition matrix would be required. Details 
of this procedure are discussed elsewhere 11811 . In the case 
of the ideal hardware error model, all erasure patterns can be 
mapped directly to encode measurements by measuring all the 
qubits of a code block that is not erasure free. In the case of 
the lossy hardware model, we take a worst case approach and 
all erasure patterns at the end of the error correction procedure 
are taken to be an encoded full erasure, so there is no need for 
the more detailed analysis - we can simply replace the cor- 
rupted block of qubits with a block in a known fixed encoded 
state. 

In the case of perfect hardware (6 = 0), the error rate recur- 
sion relation is 

eideii = 566^+4064^ + 3878e^ al - 129675ef deal + • • • , 

(8) 

which yields a threshold, for ej^j = eideai, of approximately 
Cideai = 0.115. The Markov chain describing the error recov- 
ery procedure for this model is shown in Figure|5] 

In the case of lossy hardware with perfect teleportation, as- 
suming that 6 = ei oss , the error rate recursion relation is 

e[2 = 1050ef oss +33173 e L-46242 ei 5 oss -6861701 e f oss +--- , 

(9) 

which yields a threshold, for eJJ^ = 7 eioss > of approximately 
eioss = 0.0178. This threshold is only valid if it is identical to 
or smaller than the encoded measurement threshold. Because 




FIG. 5: The Markov chain describing the error recovery for the ideal 
hardware error model. The states are labeled by the weight of the 
correctable erasures, and the 'fail' corresponds to all patterns that are 
not correctable according to our procedure. Probabilites are omitted 
for readability. 




FIG. 6: The Markov chain describing the error recovery for the lossy 
hardware error model. The states are labeled [m, n] according to the 
number of full erasures m and the number of Z erasures n, and the 
'fail' state corresponds to all patterns that are not correctable accord- 
ing to our procedure. Probabilites are omitted for readability. 



of the structure of CSS codes, encoded basis states correspond 
to superpositions of elements of a coset of a linear classical 
code. In the case of the [[7, 1 , 3]] quantum code, the linear 
classical code is a [7,4,3] code. Measurement failures can 
then be seen as classical erasures on this [7,4,3] code, and 
ignoring the correctable classical erasure patterns of weight 
three and higher, the encoded failure rate for measurements is 
given by 

§ m = W 7 Vn-5) 7 -\ do) 

i=3 

which yields the benign error threshold 6 = 0.25, validating 
the calculated threshold value ei oss = 0.0178 under the as- 
sumption 6 = 6i oss . This is a worst case assumption because 
the family of teleportation protocols described in Ref. Q| use 
an increasing number of detectors to increase the probability 
of success, and for the smallest such protocol the probability 
of photon loss is the same as the probability of photon loss for 
a single qubit measurement. The Markov chain for the lossy 
hardware error model is depicted in Figure[6] 
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VIII. CONCLUSIONS 

Using the error correction techniques outline here, the 
error threshold for Clifford gates is found to be at least 
0.0178 < e < 0.115 (since e loss < e < e idea i), where 
e is the probability that some type of erasure is introduced 
due to photon loss at the detectors or due to a teleportation 
failure. 

The threshold values calculated here can be improve d by 
using optimized stabilizer measurement techniques I22l 12311 . 
or by merging stabilizer measurement and error correction 
steps more aggressively I^. l24ll . Figure|4]is a straightforward 
generalization of a technique used for a two qubit Z measure- 
ment error correcting code |2], but Knill showed that by merg- 
ing all stabilizer measurements with error correction steps in 
modified teleportation protocols, significantly higher thresh- 
olds can be obtained 12411 . 

The Markov chain description of the error correction pro- 
cedure, discussed in more detail elsewhere fl8il . can be used 



with any of these techniques. This systematic approach to 
the calculation of the encoded error rates is particularly use- 
ful for practical applications of concatenated codes, since it is 
able to give the probability distribution of the encoded errors. 
The calculation to check the dependency of the distribution 
on parameters such as the number of error correction attempts 
or the number of concatenation levels is straightforward, and 
one could determine easily how many correction attempts or 
concatenation levels are necessary to obtain some target error 
rate. 
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